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Abstract 

Bose-Einstein correlations in W-pair production at LEP are investigated in a 
data sample of 629 pb -1 collected by the L3 detector at yfs = 189-209 GeV. Bose- 
Einstein correlations between pions within a W decay are observed and found to 
be in good agreement with those in light-quark Z decay. No evidence is found for 
Bose-Einstein correlations between hadrons coming from different W's in the same 
event. 
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Introduction 



In hadronic Z decay, Bose- Einstein correlations (BEC) are observed as an enhanced production 
of identical bosons at small four- momentum difference PUI2|- BEC are also expected within 
hadronic W decay (intra- W BEC). At LEP energies, in fully-hadronic W + W~ events (qqqq) 
the W decay products overlap in space-time. Therefore, it is also natural to expect BEC 
between identical bosons originating from different W's (inter-W BEC). A comparison of BEC 
in fully-hadronic W + W~ events with those in semi-hadronic W + W~ events (qq^), serves as 
a probe to study inter-W BEC. Together with colour reconnection inter-W BEC form a 

potential bias in the determination of the W mass at LEP. 

In this Letter we search for evidence of inter-W BEC. We examine all BEC, intra- W as 
well as inter-W, and make a comparison with the BEC observed in hadronic Z decay with and 
without the contribution of Z — > bb decay. This Letter updates the analysis of Reference U\ 
using about 3.5 times more data. 

Data and Monte Carlo 

The data used in this analysis were collected by the L3 detector jHj at y/s = 189 — 209 GeV 
and correspond to a total integrated luminosity of 629 pb _1 . These data are grouped in seven 
energy bins listed in Table [TJ 

Fully-hadronic and semi-hadronic W + W~ events are selected with criteria similar to those 
described in Reference El An additional requirement for the fully-hadronic channel is a cut on 
the neural network output [§] to further separate the signal from the dominant e + e~ — > qq(7) 
background. In total, about 3,800 semi-hadronic and 5,100 fully-hadronic events are selected. 

The event generator KORALW with the BEC algorithm BE 32 [3] is used to simulate the 
signal process. The values of the BE32 parameters are found by tuning the Monte Carlo (MC) 
to Z-decay data depleted in b-quark events. Both the BEC and the fragmentation parameters 
are tuned simultaneously This set is used in this analysis. Systematic studies are made using 
an alternative set of parameter values^ obtained by tuning to Z-decay data of all flavours and 
used in Reference HI 

The background processes e + e~ — > qq(7), e + e~ — > ZZ and e + e _ — ► Ze + e _ are generated 
using PYTHIA ^T]. For the qq(7) channel KK2f ^2] is also used. BEC are included in both 
programs. The generated events are passed through the L3 detector simulation program |13j . 
reconstructed and subjected to the same selection criteria as the data. 

The selection efficiencies of the W + W~ channels qqez/, qq/U^, qqrz^ and qqqq are found 
to be 83%, 75%, 50% and 86%, respectively. The purities of these channels are around 95%, 
95%, 85% and 80%, respectively, varying a few percent between the different energy bins. The 
selection efficiencies of fully-hadronic events changes by less than 0.5% when BEC (intra-W, or 
both intra-W and inter-W) are not simulated. 

A sample of hadronic Z-decay events is selected from the data collected by the L3 detector 
at y/s = 91.2 GeV, corresponding to a total integrated luminosity of 48.1 pb _1 . The event 
selection, which is similar to that presented in Reference El results in about one million hadronic 
Z-decay events with a purity greater than 99%. Since b quarks are greatly suppressed in W 
decays, and to compare with the W + W~ data, a light-quark (u, d, s or c) enhanced Z-decay 

^ The BEC parameter values PARj(92)=0.90, giving the BEC strength, and parj(93)=0.43 GeV, the inverse 
effective source size, are obtained. 

2 ) The alternative set uses the BEC parameter values PARj(92)=1.68 and PARj(93)=0.38 GeV. 
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sample is selected using an anti-b tagging procedure ^1] to reduce the bb fraction from 22% 
to 3%. This results in about 180,000 Z-decay events depleted in b quarks. 

The charged pions used for the BEC study are detected as tracks in the central tracker, 
using selection criteria similar to those of Reference About 82% of the tracks selected in 
MC samples are pions. This selection yields about one million pairs of like-sign particles in the 
fully-hadronic channel and about 200,000 pairs in the semi-hadronic channel. 



Analysis Method 

The Two-Particle Bose-Einstein Correlation Function 

BEC are due to interference between identical bosons which are close to one another in phase 
space. These correlations can be studied in terms of the two-particle correlation function 

Re(Pi,P2) = — r , (1) 

Po(pi,P2) 

where p2{Pi,P2) is the two-particle number density of particles with four-momenta p\ and P2, 
and po(pi,P2) the same density in the absence of BEC. As is customary in correlation studies, 
the normalization is J J ^2(^1,^2) dpi dp2 = {n(n — Su)), where n is the number of particles 
and 812 is 1 for identical particles and otherwise. The largest BEC occur at small absolute 



four-momentum difference Q = \J—{px — P2) 2 and R2 is parametrized in this one- dimensional 
distance measure by defining 

where N ev is the number of selected events and n pa i rs the number of like-sign track pairs in the 
N ev events, retaining the normalization of p2{PitP2)- 

Apart from the quark flavour, hadronic W and Z decays are expected to be similar. Hence, 
we measure the correlation function for a light-quark Z-decay sample, i?|^ udcs , in order to 
compare it with those of the fully-hadronic W + W~ events, i?^ w , and the semi-hadronic W + W~ 
events, R% . 

If there is no inter- W interference, we can write 

pf W (pi,P2) = 2pf( Pl ,p 2 ) + 2p?(p l )p?(p a ) , (3) 

where the assumption is made that the densities for the W + and W~ bosons are the same. In 
order to transform from pi,p2 to Q, we replace the product of the single-particle densities by 
a two-particle density, p mix , yielding pf w (Q) = 2pf(Q) + 2p^{Q). We obtain p^(Q) by 
pairing particles originating from two different semi-hadronic events. By construction, particles 
in these pairs are uncorrelated. The event mixing procedure is explained in detail in ReferencesE| 
and 16 This procedure also takes into account the momentum correlation between the two 
W's. The other two densities, p^ w and p% , are measured directly in the fully-hadronic and 
the semi-hadronic events, respectively. 

For comparison to the actual W + W~ data, i? 2 Arw is estimated for the case that inter-W 
BEC are absent. In this case 

RlZintJQ) = ^ww|§ = 2pw[g] + 2pw|{g] = 1 + (l - HQ)) { R f(Q) - l) , (4) 
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where 

2o ww fO) 

F(0) = mix (5) 

is the fraction of pairs, in which the two particles come from different W's, that would occur in 
fully-hadronic events if there were no BEC. In Equation (J1J) we assume that p£^ is the same 
for a sample with and without intra-W BEC (and no inter- W BEC). From a MC study, this 
assumption is found to be valid. In the extreme case of complete overlap of the two-particle 
densities, p^ = p£^? = |p^ w , and then F(Q) = 0.5. In the other extreme, namely no overlap, 
Pmb^ = at Q = 0, so that F(Q) = there, with the necessary compensation at large Q. Then 
BEC are equally strong in fully- and semi-hadronic events. Thus, depending on the degree of 
overlap, the strength of the BEC in fully-hadronic events is reduced by up to 50% compared 
to that of semi-hadronic events. We estimate F{Q) using a MC without BEC. 

Direct Search for Inter-W BEC 

The hypothesis that the two W's decay independently can be directly tested using Equation (jHJ). 
In particular, the following test statistics are defined as the difference and the ratio of the left- 
and right-hand side of Equation (jSJ) in terms of Q 

Ap(Q) = pf^(Q) - 2pf(Q) - 2p^(Q) (6) 

and 



D(O) = ^ ^ (7) 



2 P f(Q) + 2p^(Q) 

This method gives access to inter- W correlations directly from the data, with no need of MC |15j . 

In the absence of inter-W correlations, Ap = and D = 1. To study inter-W BEC, 
deviations from these values are examined at small values of Q for like-sign particles. The 
influence of other correlations or potential bias on these quantities is studied by analysing 
unlike-sign pairs and MC events. 

The event mixing procedure could introduce artificial distortions or not fully account for 
some correlations other than BEC or some detector effects, causing a deviation of Ap from zero 
or D from unity for data as well as for a MC without inter-W BEC. These possible effects are 
reduced by using the double ratio 

D'(Q) = nr ^ (9)data , (8) 

-^IVjMC, no inter 

where D(Q)mc, nointer is derived from a MC sample without inter-W BEC. 

Results 

Measurement of the Correlation Function 

To obtain the correlation function i? 2 , Equation for the W + W~ events, background is 
subtracted by replacing p2(Q), Equation ((2)), by 



where V is the purity of the selection and nb g is the number of pairs of tracks corresponding 
to (1 — V)N ev background events. This density is further corrected for detector resolution, 
acceptance, efficiency and for particle misidentification with a multiplicative factor derived 
from MC. Since no hadrons are identified, this factor is the two-pion density found from MC 
events at generator level divided by the two-particle density found using all particles after full 
detector simulation, reconstruction and selection. For this detector correction, the no inter-W 
scenario with the BE32 algorithm is used. 

The reference distribution po is calculated from a KORALW MC sample without any BEC 
at generator level, using only pions. We use a MC distribution rather than unlike-sign particle 
pairs from data, as these are strongly affected by the presence of dynamical correlations, such 
as resonances. 

Figure C] shows the correlation function for semi-hadronic and fully-hadronic events. Results 
for different energy bins are compatible with each other, and therefore are combined. 

The correlation function of the light-quark Z-decay data sample, i?|^ udcs , is determined 
similarly to and R^ w . Figure compares i?|^ udcs with R^ and, as expected, good 
agreement is observed. If b-quark decays of the Z are not removed from the sample, a depletion 
of the correlation function is observed at low Q, and a clear discrepancy exists with the W 
data, as also shown in Figure UK- 

The BEC enhancement at low Q values is fitted from Q = to 1.4 GeV with the first-order 
Edgeworth expansion of the Gaussian [H1I21IIE1, 



MQ) = 7(1 + sq) 



1 + Aexp (-R 2 Q 2 ) (l + ^H 3 (V2RQ 



(10) 

using the full covariance matrix [TB], which takes into account the effect of bin-to-bin correla- 
tions. The parameter 7 is an overall normalization factor, the term (1 + 8Q) takes into account 
possible long-range momentum correlations, A measures the strength of the BEC, R is related 
to the effective source size in space-time, k measures the deviation from the Gaussian and 
H 3 (x) = x 3 — 3x is the third Hermite polynomial. The fit results are given in Table El and 
shown in Figure d 

BEC are observed (A > 0) in both fully-hadronic and semi-hadronic events. The value of 
A is higher for the semi-hadronic than for the fully-hadronic channel, which could be due to a 
partial suppression of inter-W BEC [15J combined with incomplete overlap of p^ + and p^ . 
Using a MC without BEC to correct the data for detector effects, reduces A in both channels 
by about 30%, but the difference in the A values remains. 

Different sources of systematic uncertainties on the fit parameters are studied. Track and 
event selections are varied, the background fractions are varied by ±20%, covering biases in the 
Q distributions of background events, and different MC's, using different sets of MC parameter 
values, are used to generate the background events. Moreover, the stability of the results 
is tested by changing the fully-hadronic event selection such that the background fraction is 
reduced by more than a factor two. 

The influence of the choice of the MC used for the reference sample and for the detector 
correction is also taken into account. For the detector correction a MC with intra-W BEC is 
used, since it agrees better with the Q distribution of the raw data than a MC without BEC [7j. 
Both sets of parameter values in BE 32 are used and for the fully-hadronic channel both the 
no inter-W and inter-W scenarios are considered. For the inter-W scenario, the systematic 
uncertainty is estimated using the fraction of inter-W BEC consistent with the measurements 
of Ap, D and D' discussed below. Finally, a systematic uncertainty is estimated by varying 
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the fit range by ±320 MeV. The contributions to the systematic uncertainty of A are shown in 
Table El 

The correlation function for fully-hadronic events expected for the case of no inter-W 
BEC, -R^rm mt er ; ^ s estimated from Equation Q. In this calculation, the value of F(Q), from 
Equation (jSj) using KORALW without BEC, is around 0.2 at low Q and increases to 0.6 for 
Q > 3 GeV. Using these values of F(Q) and Rf = i?|^ udcs in Equation (jlj), we obtain the 
no inter-W prediction shown in Figure as the full histogram. Good agreement is observed 
between this distribution and the fully-hadronic W + W~ data, thus indicating no or only weak 
inter-W BEC. 

Measurement of Ap, D and D f , Direct Search for Inter-W BEC 

Figure El shows the distribution of Ap, Equation ©, for like-sign, (±,±), and for unlike- 
sign, (+,—), particle pairs. Figure El shows the distributions of D and D', Equations (J7J) 
and (JHJ), for like-sign and unlike-sign pairs. For the double ratio D' we use the no inter-W 
scenario of KORALW as the reference sample. The distributions of Ap, D and D' are not 
corrected for detector effects, but background is estimated from MC and subtracted according 
to Equation 0, from and p^ w . Also shown in Figures El and El are the predictions of 
KORALW after full detector simulation, reconstruction and selection. Both the inter-W and 
no inter-W scenarios are shown. 

The inter-W scenario shows an enhancement at small values of Q in the Ap, D and D' 
distributions for like-sign pairs. We also observe a small enhancement for unlike-sign pairs due 
to the MC implementation of BEC, which shifts the vector momentum of particles, affecting 
both the like-sign and unlike-sign particle spectra. The no inter-W scenario describes the 
Ap(±,±), £)(±,±) and D'(±,±) distributions, while the inter-W scenario is disfavoured. 

For quantitative comparisons, the integral 



is computed. Also, the D'(±, ±) distribution is fitted from Q = to 1.4 GeV, using the full 
covariance matrix, with the parametrization 



where 5, A and k are the fit parameters. Both J(±, ±) and A measure the strength of inter-W 



The systematic uncertainties on J(±, ±) and on A are listed in Tables 0] and El respectively. 
In addition to the track and event selections, the amount of background is also varied and differ- 
ent MC's, using both sets of MC parameter values, are used to generate the background events. 
Furthermore, contributions to the systematic uncertainty on A are obtained by varying the 
choice of MC for the reference sample in D' using PYTHIA and KORALW, both with no BEC 
at all and with only intra-W BEC. MC's without BEC are used to estimate the effect of residual 
intra-W BEC. The effect of various models for colour reconnection is included. Changes in 
the fit range (±400 MeV), in the bin size (from 40 to 80 MeV) and in the parametrization of 
Equation (fT2"j) (removing (1 + SQ) from the fit) also contribute to the systematic uncertainty 
on A. 

3 ^The so-called SKI (with reconnection probability of about 30%), SKII and SKIP ,6. models, as implemented 
in PYTHIA, are used. 




(11) 



D\Q) = (1 + 5Q)(1 + Aexp(-k 2 Q 2 )) 



(12) 



BEC. 
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In the mixing procedure, a semi-hadronic event is allowed to be combined with all possible 
other semi-hadronic events. To be sure that this does not introduce a bias, the analysis is 
repeated for a mixed sample where every semi-hadronic event is used at most once. The 
influence of the mixing procedure is also studied by not only combining oppositely charged 
W's, but also like-sign W's. The influence of an extra momentum [7|, used in the event mixing, 
is also included as a systematic effect. The effect of these three changes in the mixing procedure 
is also given in Tables |U and El Moreover, the analysis is repeated removing the cut on the 
neural network output for the mixed events. Furthermore, the effect of uncertainties in the 
energy calibration of the calorimeters is studied. Finally, the influence of incorrect assignment 
of tracks to r or qq systems in the qqru channel is investigated. 

The value of J(±, ±) is computed using the full covariance matrix, taking Q max = 0.68 GeV, 
the value where the two MC scenarios have converged to less than one standard deviation. The 
results for each centre-of-mass energy, displayed in Figure EK, are consistent with each other. 
Combining all J(±, ±) values results in 

J(±,±) = 0.03 ±0.33 ±0.15 , 

where the first uncertainty is statistical and the second systematic. Using KORALW with the 
inter-W scenario gives J(±, ±) = 1.38 ± 0.10, where the uncertainty is statistical only. In 
Figure E^i this value is shown as a vertical band. It disagrees with the value of the data by 3.6 
standard deviations. For unlike-sign pairs we obtain J(±, — ) = 0.01 ± 0.36 ± 0.16, consistent 
with zero. 

The value of the fit parameter A, Equation ()12|). is shown in Figure |3Jd for each energy bin. 
Combining all A values results in 

A = 0.008 ±0.018 ±0.012 , 

where the first uncertainty is statistical and the second systematic. The value of k is found to 
be 0.4 ± 0.4 ± 0.3 fm and the correlation coefficient between A and k is 0.45. A similar fit is 
performed for the KORALW MC sample with inter-W BEC, resulting in A = 0.098 ± 0.008, 
where the uncertainty is statistical only. In Figure EJo this value is shown as a vertical band. 
It disagrees with the value of the data by 3.8 standard deviations. Using the alternative set of 
MC parameters results in J(±, ±) = 1.78±0.10 and A = 0.126±0.008, where the uncertainties 
are statistical only. 

To summarize, an excess at small values of Q in the distributions of Ap(±,±), D(±,±) 
and -D'(±, ±) is expected from inter-W BEC, but none is seen. These distributions agree well 
with KORALW using BE32 without inter-W BEC, but not when inter-W BEC are included. 
We thus find no evidence for BEC between identical pions originating from different Ws. 
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(V5> (GeV) 


188.6 


191.6 


195.5 


199.5 


201.8 


204.8 


206.6 


Integrated luminosity (pb 1 ) 


176.8 


29.8 


84.1 


83.3 


37.2 


79.0 


139.1 



Table 1: Average centre-of-mass energies and corresponding integrated luminosities. 



Parameter 


Fully-hadronic 


Semi-hadronic 


7 
A 

R (fm) 
5 (GeV- 1 ) 

K 

X 2 /dof 


0.92 ±0.02 ±0.02 
0.74 ±0.04 ±0.06 
0.73 ±0.03 ±0.06 
0.00 ±0.02 ±0.05 
0.80 ±0.09 ±0.10 
27.2/30 


0.91 ±0.02 ±0.02 
0.95 ±0.08 ±0.05 
0.71 ±0.04 ±0.07 
0.03 ±0.02 ±0.05 
0.61 ±0.14 ±0.11 
35.4/30 



Table 2: Values of the fit parameters 7, A, R, 5 and k for the fully-hadronic and semi-hadronic 
events. The first uncertainty is statistical, the second systematic. The \ 2 an d number of degrees 
of freedom (dof) are also given. 



Source 


Fully-hadronic 


Semi-hadronic 


Track selection 


0.021 


0.016 


Event selection 


0.013 


0.014 


Background contribution 


0.018 


0.021 


Alternative MC as reference 


0.013 


0.013 


Alternative MC for correction 


0.022 


0.024 


Inter-W BE in MC correction 


0.032 




Fit range 


0.020 


0.018 


Total 


0.06 


0.05 



Table 3: Contributions to the systematic uncertainty on the A parameter for the fully-hadronic 
and semi-hadronic events. 



Source 


Contribution 


Track selection 


0.084 


Event selection 


0.068 


Background contribution 


0.055 


Mixing procedure 


0.065 


Neural network cut 


0.038 


Energy calibration 


0.024 


Track misassignment in qqru channel 


0.038 


Total 


0.15 



Table 4: Contributions to the systematic uncertainty of J(±, ±). 
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Source 


Contribution 


Track selection 


0.0029 


Event selection 


0.0049 


Background contribution 


0.0042 


Alternative MC as a reference 


0.0060 


Colour reconnection 


0.0026 


Fit range 


0.0018 


Rebinning 


0.0020 


Fit parametrization 


0.0017 


Mixing procedure 


0.0044 


Neural network cut 


0.0033 


Energy calibration 


0.0017 


Track misassignment in qqru channel 


0.0022 


Total 


0.012 



Table 5: Contributions to the systematic uncertainty of the A parameter. 
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Figure 1: The correlation function R 2 at \fs = 189-209 GeV, for (a) semi-hadronic events and 
(b) fully-hadronic events. The solid curves are the results from the fits of Equation (fTU|) to the 
W + W~ data. In (a) the full histogram is for the light-quark Z-decay data sample and the dashed 
histogram is for the data sample containing all hadronic Z decays. In (b) the full histogram 
gives the expectation when inter-W BEC are absent, Equation (JH), using — /?2^ u< ^ cs • 
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Figure 2: Distributions for uncorrected data at ^/s = 189-209 GeV of (a) Ap(±, ±) and (b) 
Ap(+, — ). Also shown are the MC predictions of KORALW with and without inter-W BEC. 
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Figure 3: Distributions for uncorrected data at y/s = 189-209 GeV of (a) D(±, ±), (b) D(+, — ), 
(c) D'(±,±) and (d) £>'(+,-). Also shown are the MC predictions of KORALW with and 
without inter-W BEC. 
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Figure 4: Values of (a) the integral J(±, ±) and (b) the A parameter, at different centre-of- 
mass energies and their average. The uncertainties are statistical only. The wide bands show 
the average value of the data including the systematic uncertainties. Also shown are the MC 
predictions of KORALW with inter-W BEC. 
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